Diffraction of a superfiuid Fermi gas by an atomic grating 
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An atomic grating generated by a pulsed standing wave laser field is proposed to manipulate the 
superfiuid state in a quantum degenerate gas of fermionic atoms. We show that in the presence of 
atomic Cooper pairs, the density oscillations of the gas caused by the atomic grating exhibit a much 
longer coherence time than that in the normal Fermi gas. Our result indicates that the technique 
of a pulsed atomic grating can be a potential candidate to detect the atomic superfiuid state in a 
quantum degenerate Fermi gas. 
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Since the observation of Bosc-Einstcin condensation in 
trapped atomic gases jjj there has been increasing inter- 
est in the possibility of observing the Bardeen-Cooper- 
Schrieffcr (BCS) phase transition to the superfiuid state 
Q in dilute fermionic akali gases ^ q|. Currently, 

experimental efforts in cooling of fermionic atoms of 6 Li 
Jjjand 40 K || to the quantum degenerate regime have 
made significant progress, reaching temperatures as low 
as 0.21V where Tp is the Fermi temperature. 

At the ultracold temperatures achieved in these exper- 
iments, p-wave collisions between atoms are highly sup- 
pressed and s-wave collisions between atoms in the same 
internal state are forbidden by the Pauli exclusion princi- 
ple. As a result, the most likely possibility for the forma- 
tion of Cooper pairs is an attractive s-wave interaction 
between atoms in different hyperfine states. Fortunately, 
6 Li and 40 K appear to be very promising candidates. 6 Li 
possesses an anomolously large and negative s-wave scat- 
tering length H while for 40 K, a Feshbach resonance ex- 
ists for two of the hyperfine states which can be used to 
create the required large attractive interaction H . There 
have been several recent proposals involving the molec- 
ular state formed by either a Feshbach resonance [Q or 
photoassociation |^| which would indicate that transition 
temperatures > O.lTp might be possible, thereby placing 
the BCS transition within reach of current experiments. 

Unlike Bose-Einstein condensation, the BCS transition 
is not characterized by a significant change in either the 
density or momentum distribution of the gas. Hence, how 
to identify the BCS transition and to study the proper- 
ties of the superfiuid state in these gases remain an open 
question. Recently, there have been a number of propos- 
als to detect the BCS state including off-resonant light 
scattering Jp |, measuring the frequency of low energy ex- 
citations KG |, and the line shift in the absorption spec- 
trum of a resonant laser [ p~T| . Since the Cooper pairing 
of the fermions is represented by a phase coherent su- 
perposition of pairs of single particle states in the BCS 
wave function, processes which involve the scattering of 
particles can exhibit quantum interferences that are not 
present in the normal state 12 . However, the previous 
proposals || |l0|, |ll[ are only sensitive to the existence 
of a superfiuid order parameter and/or the spectrum of 
quasiparticle excitations of the gas and do not directly 



probe the quantum coherence of the BCS state. In this 
letter, we study the possibility of using atom optics tech- 
niques to directly study the coherence presence in a gas of 
fermionic atoms that have undergone the BCS transition. 



In our proposal, a pulsed atomic grating is generated 
by a far detuned standing wave optical field of frequency 
u>l and applied to the Fermi gas, similarly to the recent 
experiment with a Bosc-Einstcin condensate by Deng et 
al. ]l3|. The grating results in the scattering of atoms 
from states of momentum Kk to S(k + iq) where ±q/2 
are the wave vectors of the travelling waves forming the 
grating and I is an integer. For |q| <C 2kp, there is 
a rapid decay of the density oscillations in the normal 
state due to a dephasing caused by the range of energies 
near the Fermi surface that are excited by the grating. In 
contrast, for the BCS state, even though this dephasing 
also initially reduces the amplitude of the density oscil- 
lations, small amplitude oscillations nevertheless persist 
indefinitely. They result from the constructive quantum 
interference for the creation of pairs of quasiparticles by 
the grating and the large density of states for low energy 
quasiparticles. The detection of such persistent density 
oscillations would provide a clear indication of the BCS 
state. 



We consider a box of volume V containing A" fermionic 
atoms with density n — N/V . The gas is assumed 
to be at zero temperature (T = 0). The atoms pos- 
sess two degenerate ground states, |f) and ||), that have 
equal population and equal chemical potentials, fi. For 
convenience, we refer to the two states as spin states. 
The fermion annihilation (creation) operators Ck a (cL a ) 
destroy (create) an atom in a state with momentum 
fik and spin \a) and obey the anti-commutation rela- 
tions jcko-,4'0-'} = <*k,k'<W' and {cko-jCW} = 0. The 

Hamiltonian for the system is given by H = Ho + Hi 
where Hj represents the interaction of the atoms with 
the optical potential formed by the laser and Hq is the 
free Hamiltonian for the fermions. It includes attractive 
s-wave interactions between atoms in opposite spin states 
with scattering length a < 0. By introducing the order 
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parameter for the BCS state, 

. ATth 2 \a\ „ 

~^r^ (c - kiCkT) ' 
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-ffo rnay be diagonalized by a canonical transformation 

IBS, 
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other, thereby resulting in a momentum transfer to the 
atom of ±q. 

We wish to examine how the atomic density, 
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q' 



evolves in time following the application of the grating. 
We work in the Heisenberg representation with the time 
dependence for the operators given by 



where Eq is the ground state energy and a ka (d k(T ) are 
quasiparticle creation (annihilation) operators given by, 

d kT = cos(6» k /2) c kT - sin(6» k /2) cL ki , (3a) 
d f _ ki = cos(0 k /2)c! Ui +sin(0 k /2)c kT . (3b) 



The energies of the quasiparticles are fii/ k = + A 2 
while tan0 k = A/£ k . Here, £ k = h 2 k 2 /2m — £p is the 
kinetic energy measured relative to the Fermi energy, 
ep = ti 2 k F /2m = n + 2Trh 2 \a\n/m. The ground state 
of Ho, |$) , is the vacuum state for the quasiparticles, 
dko- 1 3?) =0. For A = one recovers the normal state of 
a Fermi gas subject to a Hartree-Fock mean-field energy 
of — 2irh 2 \a\n/m for each spin component. 

At time t — a standing wave laser field is applied 
to the gas for a duration r, coupling the atoms to an 
excited state of frequency u> a with a Rabi frequency Q. 
For detunings 6 = ojl — to a large compared to the excited 
state linewidth and f2, the excited state may be adiabat- 
ically eliminated resulting in an optical dipole potential 
Uq cos 2 (q • x/2) where Uq = HQ 2 / 48. In the Raman-Nath 
regime, the energy of the quasiparticles excited by the po- 
tential may be neglected in comparison to the interaction 
energy, i.e., U > ^k F +k max where hk max = U rq is the 
maximum momentum imparted by the potential in the 
time r Q|. For the condition q <C 2Uf considered in this 



paper, this requires r <C 



l-(A/(7 ) 2 /(hk F q/m) 



In this limit, the potential acts as an impulsive phase 
grating with the second quantized form, 



n I = -h\5(t) (p q + /S_ q ) , 

where A = —Ttt 2 /168 and 



pi = 4 



k+qo- Ckcr 



(4) 
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is the q th Fourier component of the atomic density. Phys- 
ically, Eq. (|4|) may be interpreted as the result of an atom 
absorbing a photon from one of the travelling waves fol- 
lowed by the stimulated emission of a photon into the 



c kCT (t) 



-i\(p n +p-q) e iH a t c ko . e - i - f/ °* e iA (Pq+P-q). (7) 



A direct application of Eq. (0) to calculate c kcr (t) for 
t > would, in general, lead to incorrect results be- 
cause Eq. (g) depends on the mean-field A, which is 
modified by the scattering off the optical potential and 
must be calculated in a self-consistent manner. For 
t > 0, one must use the new value of the order pa- 
rameter to perform a canonical transformation on Hq. 
Let us define the pair expectation value at the time 
t = + immediately after the application of the grating 
c_ k < i (0 + )c kT (0 + )). Using the result [fijl, 



as Fitl 



iV s (2A)c k _ s 



(8) 



where J s is a Bessel function of integer order s, we have 
the explicit form for F^l, 



1 °° 

*i£k=9 * rl ^(2A)J„- s (2A)sin# : 



k— sq^k— k' ,nq- 



(9) 



The mean-field following the grating, A^ + ) = 

47 y"J, a S k k' F^~-L e ~ l< - k ~ k '' x , is now spatially modulated 
at all harmonics of q. Since the probability that a pair 
state (k T , — k J.) is scattered by the grating into the 
pair state (k + nq | , — k - mq j.) is J 2 (2A) J^(2A) ~ 
A 2 (l™l+l m D, for A < 1 the mean-field experienced by each 
atom will not differ significantly from what it was for t < 
0. Therefore, we restrict ourselves to A < 1 and use the 
zeroth order approximation A^ + ) « A. In the examples 
presented below, we have numerically calculated F^^ + 
and confirmed that for n ^ 0, F^ +n jF^ < 0.1. 

Using Eq. (^), the scattering of the atoms at t = 
may be treated to all orders in A. After some algebra 
one obtains an expression for the density, 
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p(x,i) = Po + ^Y. £ J s (2A)J n+s (2A){/ s , n (t)cos[(q-x + 7r/2)n]-ff s>n (t)siii[(q-x + 7r/2)n]}, (10) 

n=l s— — oo 

i 



i 

fa,n(t) = ^ ^ Sin jjflk-sq I COS I^t COS Vu+ngt Sm jjflk-sq + SU1 V*t SU1 Z/k +wq tyl s , w (k) 1 , (11) 

&>,«(*) = — sin ^^k-sg {sin v\t cos iVk+ng*-B s ,o(k) - cos z/kt sin ^ +wq t-B S|n (k)} , (12) 

k 
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FIG. 1: Density oscillation of the Fermi gas after the grat- 
ing, (a), (b) The spatio-temporal evolution of the atomic 
density; (c) density oscillations at x = 0. The curves for 
A = 0.1 and 0.2 have been shifted down by 0.15 and 0.3, re- 
spectively. The units for time is the so-called Talbot timejl^] 
tt = 2-ir/[hq 2 /(2m)], and the units for energy and density are 
ef and /cf, respectively. Here A = 0.2, q — 0.3&;f- For (a) 
and (b), the brighter regions represent higher density. 



where p is the spatially uniform background density, 

A s>n (k) =sin [i0k-s q - (-i) s flk - (-iy +n e k+nq \ + 

(-1)" sin 0k sin6> k +riq/ sin |0k-sqi and 
B.,„(k) =sin [i0k- S q - (-ir +n ^k+nq]- For A = 0, we 
recover the results of Ref . Ea| . 

Since the formation of Cooper pairs by the atoms 
only affects those states lying in a shell of thickness 
8k fts mA/h 2 kp around kp, q must be chosen such that 
q ~ 5k in order to be sensitive to the quantum coher- 
ence near the Fermi surface. This condition can be reex- 
pressed as q/kp ~ A/ep. For 6 Li or 40 K, realistic val- 
ues of the zero temperature gap are A/ep ~ 0.04 — 0.1 



I 

H U , although recent models that account for the forma- 
tion of a molecular condensate formed using a Feshbach 
resonance or photoassociation predict A/ep ~ 0.1 — 0.5 
H |(| . Fig. |l| shows the density evolution for various A, 
with q along the x-axis and the k-space of the fermions 
taken to be one-dimensional for simplicity. We observe 
that the initial oscillations in the density die out in a 
characteristic time td for both the normal state and su- 
perfiuid state. This decay is attributable to the single 
particle states near the Fermi surface, with the range 
of energies SE that are excited by the grating. The 
decay time is therefore roughly given by td ~ H/5E. 
For A < 1, q is the dominant spatial harmonic in the 
density. Consequently, for the normal state, only those 
states lying within a distance q from the Fermi surface 
can be excited by the grating due to Pauli blocking so 
that 8E « — £k F -q| and td ~ m/hqkp. For the 
BCS state, those initial states lying within the shell 8k 
of kp give the greatest contribution when excited by the 
grating and therefore, SE w h\vk. p — fk F -5k| ~ 2A and 
td ~ h/2A. The ratio of dephasing times for the normal 
state and the superfiuid state is {A/ep) (q/kp)^ 1 ~ 1. 
This explains why the decay times in Fig. ffl(c) appear to 
be the same. 

The major difference between the superfiuid state and 
the normal state is the persistence of almost undamped 
density oscillations for t > td in the former case. From 
Fig. |l], we see that for the normal state, the density os- 
cillations collapse at times equal to lrp/2 (I is an inte- 
ger) with incomplete revivals between collapses. This is 
due to the discontinuity in the momentum distribution of 
the normal state at kp, which leads to a dephasing term 
oscillating in time. For finite temperatures, the momen- 
tum distribution is a smooth function of k in the normal 
state and hence revivals will not occur. Such collapses 
and revivals are absent for the BCS state even at T = 
since the momentum distribution in the ground state is 
a smooth function of k for all k. The origin of the un- 
damped oscillations may be understood by evaluating the 
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FIG. 2: Oscillation frequency of the density modulation as a 
function of A. the solid line is the analytical result (see text) 
and the stars are numerical results obtained from Eq.([to|). 
The units for frequency and energy are £f/S and sf, respec- 
tively. Here A = 0.2, q = 0.2k F . 



density to lowest order in perturbation theory, 

p(x, t) = p + yYl sin2 \ ^ k +i + 0k ) 

sin [(^ k +q + fk) t] cos(q • x). (13) 

The physical effect of the grating is to excite pairs of 
quasiparticles with energies ftt'k+q and hv^. However, 
when Eq. (0) is expressed in terms of the quasiparticle 
operators, one finds that both ck+qtCVf and clk^c-k-qx 
contain the term ak+qT^ki with the respective ampli- 
tudes cos (#k+ q /2) sin (0^/2) and cos (6k/ 2) sin (#k+ q /2) . 
As a result, the transition amplitudes for the two pro- 
cesses in which an atom with spin ||) is scattered by 
the grating from k to k + q and an atom with spin J.) is 
scattered by the grating from — k — q to — k add construc- 
tively to give a total probability of creating a quasipar- 
ticle pair proportional to sin 2 [(#k+ q + #k) /2] • The con- 
structive interference reaches its maximum when k + q| 
and |k| are within Sk of kp so that #k+q ~ $k ~ 7r/2. In 
contrast, for the normal state, the quasiparticle pair cor- 
responds to an atom excited above the Fermi surface and 
a hole inside the Fermi sea and the terms c k + q fCkt and 

c^ k |C_k-qX are independent of each other. Alternately, 
one may note that if |k| and |k + q| lie within Sk of /cp, 
then the pair states (k f , k j) and (k + q f , k — q J.) 
have formed Cooper pairs in the supcrfluid state. Conse- 
quently the terms c k+q f Ckf and c^kjC-k-qi are no longer 
independent and exhibit a quantum interference due to 
the coherence between the single particle states in each 
of the Cooper pairs. 

In order to account for the degeneracy of the quasi- 



particle states with energy e = hv^, the summation over 
k may be replaced with an integration over the quasi- 
particle energies. The density of states for quasiparti- 
cles with energies in the interval [e, e + de) is N s (e)de 

- (e 2 - A 2 y 1/2 N n (£)0(e-A)de where 2V n (£) is the den- 
sity of states in the normal state|l2|]. Since N s (e) be- 
comes very large as e — > A jL(|, the states which give 
the greatest contribution to Eq. ( p"3] ) are those which 
lie on the Fermi surface. This corresponds to the terms 
in Eq. (|l3| ) with frequency ^k F +q + ^k F , which may be 
separated out from the sum. The remaining states (with 
k| ^ kp ) in the summation over k, which get out of 
phase with each other in a time on the order of td, lead 
to the collapse of the initial large amplitude oscillations. 

Consequently, the undamped oscillations are a result 
of the combined effect of the constructive interference for 
the creation of quasiparticles near the Fermi surface and 
the singular nature of the density of states for quasipar- 
ticles with momentum near the Fermi surface. A similar 
effect is observed in the relaxation rate of polarized nu- 
clear spins in superconductors, which shows a sharp rise 
for temperatures just below T c JT^ . 

Figure || shows a comparison between fi(^k F + q + 



»k F ) « A+ 



■Aep(q/k F ) 



for q/kp <C 1 and the 



oscillation frequency obtained from the numerical eval- 
uation of Eq. ( |To| ) for t ^ td- Although there is some 
discrepancy, the analytical result, ^k F +q + ^k F , does ap- 
pear to predict the correct dependence of the oscillation 
frequency on A. 

To put these results in perspective, we consider 40 K 
with a density of n — 2 x 10 14 cm~ 3 and a — — 53nm [0, 
which gives kp = (3TT 2 n) 1 ^ 3 = 1.8 x 10 7 m -1 . The con- 
dition |q| <C 2kp then corresponds to a laser wavelength 
satisfying Xl ~S> 350nm. Since the spatial extent of the 
trapped Fermi gases is typically about 0.5mm the gas 
can be treated as locally homogenous for a laser in the 
near infrared. Using A = 0.49ef exp(— 7r/2fcp|a|) p7| , 
we have A = 0.095£f. This gives a period for the un- 
damped oscillations, Ta ~ h/2A, of around 0.1ms which 
is much less than the typical life time of the trapped gas. 

In conclusion, we have shown that the manipulation 
of a degenerate Fermi gas using atom optical techniques 
will be a sensitive probe of the quantum coherence of the 
BCS state. 
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